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^ Abstract 

A method to derive stationary solutions of the relativistic Vlasov-Maxwell system is explored. In 
(~| the non-relativistic case, a method using the Hermite polynomial series to describe the deviation 

from the Maxwell-Boltzmann distribution is found to be successful in deriving a few stationary 

s 

solutions including two dimensional one. Instead of the Hermite polynomial series, two special 

orthogonal polynomial series, which are appropriate to expand the deviation from the Maxwell- 

Q Jiittner distribution, are introduced in this paper. By applying this method, a new two-dimensional 

• 1— I 

^..^ equilibrium is derived, which may provide an initial setup for investigations of three-dimensional 

relativistic collisionless reconnection of magnetic fields. 



> 
00 
0\ 
OO 

o 

\6 
o 

OO 

o 

> 



PACS numbers: Valid PACS appear here 



* Electronic address: suzuki@resceu.s.u-tokyo. ac.jp 
^Also at Department of Astronomy, Graduate School of Science, University of Tokyo, Bunkyo-ku, Tokyo 

113-0033, Japan 



1 



I. INTRODUCTION 



Relativistic plasmas play important roles in high energy phenomena. In nature, gamma- 
ray burst afterglow is considered to be a phenomenon resulting from interaction of ultrarel- 
ativistic plasmas with a Lorentz factor 7 ^ 10^ — 10^ [H E]. Detailed observations of high 
energy astrophysical phenomena, such as supernova remnants, gamma-ray bursts, and active 
galactic nuclei, reveal that non-thermal components in their spectrum comprise of radiation 
from relativistic particles. In recent laser experiments, interaction of laser and plasmas is 
realized in highly relativistic regime. These phenomena often involve sufficiently rarefied 
plasma that the effect of collisions is negligible. In order to give plausible explanations for 
such phenomena, many authors have devoted their energies to investigations into properties 
of relativistic plasmas, i.e., instabilities, radiation, and stationary equilibria. 

Here I propose a new method to obtain stationary equilibrium configurations of colli- 
sionless plasmas described by the relativistic Vlasov-Maxwell system. Such solutions are of 
great interest because they give self-consistent configurations of relativistic plasmas, which 
are used as an initial setup for studies on the relativistic coUisionless reconnection of mag- 
netic fields in El E]. In some recent numerical simulations of the relativistic coUisionless 
magnetic reconnection [3, [8], relativistic extensions of the harris sheet equilibrium [9] are 
used as the initial setup. While stationary solutions of the non-relativistic Vlasov-Maxwell 
system are well explored PIinillllinillSlElIISllISlIIZlIIH], exploration of the relativistic 
Vlasov-Maxwell equilibria is not sufficient. Although Ref. [12] discussed the existence of 
such solutions and derived a special class of them, studies on the derivation of the concrete 
expressions of the relativistic Vlasov-Maxwell equilibria are still rare. 

In the previous paper j20], a novel method to construct stationary solutions of the non- 
relativistic Vlasov-Maxwell system is proposed. In the method, the key concept is to describe 
the deviation of the distribution function from the Maxwell-Boltzmann by orthogonal poly- 
nomial series. The Hermite polynomial series is chosen for the purpose, because the weight 
function is Gaussian, which is equivalent to the Maxwell-Boltzmann distribution. Ref. [20] 
reveals that this choice is appropriate to recover some well-known equilibria (the Harris 
sheet and the Bennet pinch) and derive a new equilibrium. In this paper, I will extend 
the previous method to deal with relativistic plasmas. Instead of the Maxwell-Boltzmann 
distribution, I use a relativistic extension of the Maxwell-Boltzmann distribution, i.e., the 
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Maxwell- Jiittner distribution [21] as the weight function of these polynomial series. To the 
best of my knowledge, there is no appropriate orthogonal polynomial series that describes 
the deviation from the Maxwell- Jiittner distribution as the Hermite polynomial does for 
non-relativistic plasmas. So I introduce two orthogonal polynomials in this paper. Applying 
the method, I can derive a new two-dimensional equilibrium, which is a relativistic extension 
of the equilibrium proposed in the previous paper. This equilibrium may provide an initial 
setup for numerical simulations of the relativistic coUisionless reconnection. 

This paper is organized as follows. In Sec|Tl| I describe the procedure of the method. In 



SecTII, I apply the method to derive a stationary solution that I can treat it analytically. I 



conclude this paper in Sec, IV Appendix gives some algebraic relations used in this paper. 



II. FORMULATION 



The relativistic Vlasov equation describes the kinetic evolution of the distribution 
function /j(t,x, p) of particles j (= i for ions and e for electrons) in the phase space 
{x,y, z,px,Py,Pz) and the Maxwell equations describe the evolution of the electromagnetic 
fields. Here the cartesian coordinates in the real space are (x, y, z) and the corresponding 
coordinates in the momentum space are {Px,Py,Pz)- This system describes the exact behav- 
ior of relativistic coUisionless plasmas. In this section, I derive stationary configuration of 
plasmas uniformly extending in the z direction governed by the relativistic Vlasov- Maxwell 
system under the assumptions described in the following subsection. 



A. Equations 

Because physical variables do not depend on t or z, the relativistic Vlasov equation is 
expressed as 



+ |p|2 9^ ^m]c^ + |p|2 9y ^m]c^ + |p|2 ^ ^P- 

+ ^ \ E + q\ ^ + £^ [ P- B - Pi/ bA ^ = 

C I . /m2c2 + |p|2 j 9py C \^^^2^2 + |p|2 ^ ^^2^2 + |p|2 "j dp, 

(1) 
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where fj is the distribution function for particles j with the charge qj and the mass mj, c 



is the speed of hght and |p| = a/p^ + + Pz the norm of the momentum. Ex,Ey and 
Bx,By represent the x, ?/-components of electric and magnetic fields, which are functions of 
X and y. These electromagnetic fields are written as 



^ E = - — 
dx ' ^ dy 

^^^^ 

by introducing the scalar potential (j){x,y) and the 2;-component of the vector potential 
Az{x, y). The other components are assumed to vanish. These potentials satisfy the Poisson 
equations; 

^jli + ^jli - _47rj (3b) 

where p{x, y) and jz{x, y) are the charge density and the ^-component of the electric current 
density, respectively. These are expressed in terms of fj{x,y,Px,Py,Pz) as 



P 



/ dpx / dpy / dpjj, (4a) 



OO •/ —CO J — oo 

OO /"OO /"OO 



X^^j / dpx I dpy I dpz ^=fj, (4b) 



m^c^ + IpI 



which close the system. Solutions of Equations ([T])-(|4]) give self-consistent configurations of 
collisionless plasmas. 



B. Derivation of stationary solutions 

In the following, I describe the procedure to derive solutions of Equations ([T])-(|4]). The 
non-relativistic counterpart derived in the previous paper [20] was based on the assump- 
tion that the deviation from the Maxwell-Boltzmann distribution can be expanded by the 
Hermite polynomial series. For a relativistic plasma, the distribution function in the ther- 
modynamical equilibrium is given by the Maxwell- Jiittner distribution a relativistic 
extension of the classical Maxwell-Boltzmann distribution; 
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(5) 



where ks is the Boltzmann constant. Tj and nj represents the temperature and the density 
of the particles j. I have assumed they are constant. K2{Cj) is the modified Bessel function 
of the second kind. Some dimensionless variables have been introduced as Q = rrijC^ / {k-oTj), 
Px = Px/{mjc), py = py/{mjc), and p^ = pj{mjc). 

In seeking solutions of Equations ([T])-(|4]), I assume that the distribution function takes 
the following form; 

= IE [9]Ti^^)Sr^P.) + (p.)] I exp (-^) ff\ (6) 

Ln.=o J ^ B i/ 

where 5^™° and 5',°'^'^ are orthogonal polynomial series defined in Appendix |b} Their coef- 
ficients f?!™" and 5'°'^'^ are assumed to be functions of Az. The distribution function in this 
form describes the deviation from the Maxwell- Jiittner distribution in terms of the orthogo- 
nal polynomial expansion. The spatial dependence of the distribution function is described 
through those of potentials and A^. For large these polynomials can be reduced to the 
Hermite polynomial series as 



^even ^ -^2n(v2Cp^) ^ -ff2n+l ( V 2Cp^) 

Thus this method becomes equivalent to the previous method [2^ in the non-relativistic 
limit. 

At first, substituting the expressions for fields ([2]) and the distribution function ^ into 
Equation (jl]), we obtain 



E 



J even dn°^'^ n ■ o'^^^^ rl R^^'^'^ n o^^^fl'^" 



= (8) 



dAz " dAz " rrijC^ dpz rrijC^ dpz 

Then, to obtain equations to determine the coefficients gj^^^ and g]^, both sides of this 
equation is multiplied by S'^^'^{pz) and integrated with respect to p^, py, and pz- The 



orthogonality relation of S'f^^^ (BIO), the relation (B19a), and the expression (B21) yield 



the following ordinary differential equations 



„even 
Jj,m 



oo 



J2Bf^9]f = 0. (9) 



dAz rrijC 

■' k=m 



On the other hand, multiplying both sides of Equation (|8| by S!;^^{pz)/ ^/l + Pi + Pi + Pi 
and integrating with respect to px, py, and pz, other equations for coefficients 



if-^ E ™°-0, (10) 



fe=m+l 



are obtained. Here the relations (B15) and (B19b) have been used. Equations ^ and 
(10) give the relation between gjj^^ and 9°^- However, it is not practical to directly solve 
these infinite number of equations. Some restrictions are needed in order to solve them. 



For example, the condition (^l™" 



odd 



0{n> N) or = c/|™^i = 0(n > A^), reduces 



Equations (|9j) and (10) to a finite number of ordinary differential equations and enable us 
to obtain expressions of and as functions of A^. 

To derive the source terms in the Maxwell equations, substituting Equation ^ into 
Equations ^ and using the orthogonality relations (BIO), (B15), and (B19) again, these 
terms are written in terms of the potentials as 



niQigifl exp 



"^ili odd 

-9ifl exp 



kBT^ 



+ neqeglT exp 



I ^eQ'e odd 

+ —^9ifl exp 



Qe 



Ce 



ksTe 

kjiTe 



(11a) 
(lib) 



where the expressions (B6 ) and (B9 ) have been used for evaluation of integrals. Substitution 
of these expressions into Equations ^ leads to 



2 



Qy2 



-An 

-An 



^iQigTr exp 



^il-i odd 

9lo exp 







k^iTi 
kBTi 



+ neQegTo'' exp 



I ''^eQe odd 

+ —9e,o exp 



Qe 



Ce 



ksTe 
ksTe 



(12a) 
(12b) 



Since gjj^^ and g°^ are functions of A^, these equations determine the potentials and A^ 
as functions of x and y under certain boundary conditions. 

As well as the previous method [20], there is a limit in this method. For a stationary 
equilibrium, distribution functions can take various forms as long as the pressure balance is 



achieved. However, Equations ( 12 ) look as if there is a one-to-one correspondence between 
the field configuration and the distributions g^o^^ and g^'o^- This disagreement arises from 
the assumption that the deviation of the distribution function from the Maxwell- Jiittner dis- 
tribution can be expanded by the polynomial series S*^™"^ and 3°'^'^. Therefore, the stationary 
solutions derived above cover a part of many possible equilibria. 



III. APPLICATION 



In this section, I consider an application of the stationary solutions derived above to a 
plasma in a charge neutrality comprising of electrons and ions with the same charge but 
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the opposite sign {qi = —q^ = e and Ui = Ue = no) in which the electric field strength is 
sufficiently small (0 = 0). Then, Equation (12a) implies that ions and electrons have the 
same spatial distribution; 



even even 

yi,0 ~ iJefl ■ 



(13) 



As mentioned in the previous section, to solve Equations ([9]) and (10), I truncate both of 



,odd 



the series as f?!™" = 0(n > 2) and g"^ = 0(n > 1). Then, Equations (|9|) and (|10|) reduce to 



the following three ordinary differential equations; 



dg 



even 
J,0 



od 



„odd 



TTljC 



^godd 



2q 



dA, 
It 



TTljC 



3 ^even 



(14) 



which have solutions in the form of 



.ven =^even(Q) 



dA. 



rrijC^ 



+^?fo'(o)^+4r(o), 



^odd^2^7X(0)!^ + ^7.-?(0), 



(15) 



rrijC^ 



9]T =cr(o)' 



where (7|™"(0), (7°o'^(0), and gjf^{0) are constants of integration, which can take various 
values as long as the distribution function fj is positive at arbitrary points in the phase 



space and the condition (13) is satisffed. For example, I assume that they take the following 
values; 

0, cr(o) = i' 



gtr{0) = C, ^°^^(0) 

^:r(o) = c, gtfio) 



0, cr(o) 



m 



m. 



e_ 
2 ' 



(16) 



where C is a dimensionless constant. Some algebraic manipulations results in the distribution 
functions in the form of 



fi 



no 



fe 



2nckBTiK2{Ci) 



no 

27lckBTeK2{Ce 



eA, 

mfc^ \ c 



X exp 



+ C 



0^2(0) 



Pz + 



X exp 



eA, 



mjc^ +pI+pI+pI 



(17a) 



mlc^ +pI+pI+pI 



(17b) 
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These exactly satisfy the relativistic Vlasov equation ([T]). In order that the distribution 
function of ions /j takes positive values, the conditions 

are required. 

Next, the equations governing the field configuration is deduced by substituting the so- 



lutions (15) and (16) into Equation (12b) 



d Az d Az Sire ( nivrti njUe i , „s 
+ = V- + (19) 



Q^2 Qy2 

This equation has a solution in the form of 

Az = -BoL [cos(x/L) + cos{y/L)] , (20) 
where L is a constant satisfying the following relation; 



L = — I ^1^5 (21) 

2e y 27r(rajmjCe + neinXi) ' 

This vector potential generates a sinusoidal magnetic field as 

= Bo sm{y/L) , By = -Bo sm{x/L) . (22) 



Equations (17) and (20) provide a self-consistent configuration of relativistic plasmas. It 
is found that this equilibrium is a relativistic extension of the equilibrium proposed in the 
previous paper [20] . 

Figure [T] illustrates the thus derived configuration of equilibrium. The gray scale and the 
arrows represent the density distribution g^™ of electrons or ions in arbitrary units and the 
magnetic field, respectively. As in the non-relativistic counterpart [20], we can see that the 
current filaments lie along the z-axis and generate the magnetic fields around themselves. 
Each filament is surrounded by four filaments that carry anti-parallel currents. 

Here, I focus on the momentum distributions at two points. One is the center of a filament 
(referred to as the 0-point) and the other is a middle point between the filaments carrying 
parallel currents (referred to as the X-point). The magnetic field at the X-point is sheared 
like the Harris sheet equilibrium. Figure [2] shows the momentum distributions of ions at the 
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X and 0-points. The values of parameters are as follows; Bq = O.lmeCUe/e, where Ue is the 
electron plasma frequency, and fcTj = kT^ = mic^. From these values, one obtains 



So I have assumed C = 0.0005. While a symmetric double-peak distribution appears at the 
X-point due to = 0, an asymmetric one is achieved at the 0-point. 

In contrast to ions, the momentum distribution of electrons has a symmetric double-peak 
both at the X and 0-points, because the term e^Al/mfc^{< O.QAmlculL'^ /mfc'^) is much 
smaller than C{= 0.0005). 

Finally, I make a remark on one of the relativistic effects. In comparison with the non- 
relativistic counterpart, the distribution tends to have a relatively shallow slope in the high- 
energy regime, which comes from the difference in the equilibrium distribution. While the 
Maxwell-Boltzmann distribution is proportional to exp(— p^/2mjA;BT), the Maxwell- Jiittner 
distribution is proportional to exp{—cpz/kBT) in the limit of Pz/rrijC ^ 1. In other words, 
the relativistic one has more high-energy particles than the non-relativistic does. 

IV. CONCLUSIONS 

In this paper, I have developed a novel method to construct stationary solutions of the 
relativistic Vlasov-Maxwell system. By applying the method, a new two-dimensional equi- 
librium, which is a relativistic extension of the previous work [20] , is proposed. A comparison 
of the non-relativistic and the relativistic equilibrium is done. As well as the non-relativistic 
equilibrium, a sheared magnetic field is generated in the relativistic one. On the other hand, 
they show different behaviors in the high-energy regime. It may provide an initial setup for 
investigations of the relativistic collisionless magnetic reconnection that takes into account 
three-dimensional effects. It will be intriguing to investigate the stability of the equilibrium 
derived in this paper. 
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APPENDIX A: PROPERTIES OF THE MODIFIED BESSEL FUNCTIONS OF 
SECOND KIND 

In the next section, I often use the properties of the modified Bessel functions of second 
kind. Therefore, I review the properties in this section [22] • 

The nth-order modified Bessel functions of second kind Kn{x) is defined by the following 
integral; 

/■oo 

Kn{x)= e-^^°'^^cosh(ne)rf^. (Al) 
Jo 

It is known that Kn{x) and -ft'n+i satisfy the following recurrence formula; 



After some algebraic manipulations, Kn{x) is expressed in another form as 



(A2) 



K^{x) =('^y ^^^^^f,^, r sinh^" tdt. (A3) 



.2J T{n + 1/2) jQ 

Dividing both sides of this equation by x" and differentiating with respect to x yields the 
relation 

ir„+i(x) = (I)" r(!i+^l/2) ^°°g"''°^''coshtsinh^"tcit. (A4) 
In the limit of x ^ 1, Kn{x) reduces to 

Kn{x) ^ ^e-\ (A5) 

APPENDIX B: CONSTRUCTION FO ORTHOGONAL POLYNOMIALS 

In this section, I construct two special orthogonal polynomial series used in this paper. 



1. Some useful integrals 

In constructing the polynomials, two integrals are used again and again. These integrals 
are introduced beforehand. 



10 



The first one is 

jodd _ 



OO PCX) PCX) 



— OO •/ — OO J —cx 



^2n 

Pz exp 



-CJi+Pl + Pl + Pl 



dPxdPyPz 



(Bl) 



^/l + PI + PI + pf 

where n is an integer and ( is an independent variable. Introducing the spherical coordinates 

(p, e, 0) as 

= p sin 6^ cos 0, Py = p sin 9 sin (p, pz=p cos 9, (B2) 
the integrations with respect to 9 and can be performed to obtain 



rodd 



An 



2n+ 1 



p2"+2exp 



dp 



1+p' 



Furthermore, a change of the integral variable from p to t defined as 



p = sinht, 



leads to 



rodd 



An 



2n + l Jo 



Using the relation (A3), the integral is expressed in terms of Kn+i as 

^„dd 2"+27rr(n + 1/2) 



C"+ir(i/2) 



-^n+l(C) 



The second integral is 

jevcn _ 



OO PCX PCX 



cx J —cx J — OO 



PT exp 



<xll+Pl+Pl+Pl 



dp^dpydp^ 



Again, from the relations (B2) and (B4), this integral is converted to 

An 

e-^'=°'^^*sinh2"+2^coshtcit 



2n + 1 Jo 



Using the relation (A4), the integral is expressed in terms of Kn+2 as 

2"+27rr(n + l/2) 



reven 
n 



C"+ir(i/2) 



-^n+2(C) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



2. Even orthogonal polynomials 

I construct the polynomial series S^'^^{pz)- This polynomial must satisfy two conditions; 
(1) These are even functions. (2) The set 5*^™° forms an orthogonal basis in the even function 
space, in other words, each polynomials satisfy the following orthogonality relation. 



OO POO POD 



OO ^ — OO ^ — OO 



ST''iPz)Sr''iPz) exp -CJl +PI+PI+P 



.xdpydp^ oc 5„ 



(BIO) 
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where I choose exp [— + f)i+ f)y+ p^] as the weight function and 5mn represents Kro- 
necker's delta. Then we define the polynomials S^^'^ as 



Qeven 



oeven 



P 



2n 



n-1 
A:=0 



even ceven 
nk ^fc ) 



(Bll) 



where are coefficients determined by the orthogonality relation (BIO). Multiplying 



S^^^ in (Bll) by S'^'^'^exp [— ^ ^/l + Pl + + Pl] (rn < n) and integrating it with respect 
to px, py, and pz, the following expression for A"^™ is obtained; 

.even ^ _ I-oo IZo IZo P^'^r" ^^P ["C y^l + PI + PI + Pi] 

n n i-jsT'^r^^p Hv^+Pi+Pi+Pi] dPxdPydPz • ^ ^ 

Using this equation, one can determine the nth order polynomial S"^™"^ by induction. For 
example, the denominator and numerator of A'^'q'^ become /q™° and I'^'^"^^ respectively. Thus 



yj^even^ ^even 



Pi) is expressed as 



4 even 
^10 



A'3(C) 



amy 

In such a way, the even polynomials 5"^™° are obtained from the lower order polynomials 
m <n). The expression of the first three are 



ceven 



oeven 
•^0 



p. 



^even 



Pz 



C^2(C)' 

4 3 5i^5(C)i^2(C)- ^4(0^3(0 

C m,{C)K^{C)-[K,{C)Y 



Pz 



3i^4(C) 



(B14a) 
(B14b) 

(B14c) 



3. Odd orthogonal polynomials 



The odd polynomial series S'^'^'^{pz) are supposed to satisfy the following orthogonal 
relation, 

dpxdpydpz 



00 PCX) POO 



—00 J —00 J —00 



^r(p.)5r(p.)exp 



<\/^+Pi+pi+pi 



^+pi+pi+p'^ 



oc 6„ 



(B15) 



Note that the weight function differs from that of the even polynomial series. Then I deffne 
the polynomials 5',°'^'^ as 



n-l 



oodd _ - Qodd _ "2n+l , \ ^ /loddQi 

^0 — Pz, '-'n — Pz "T / ^nk ^/ 



odd 
k ) 



(B16) 



k=0 
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where are coefficients determined by the same procedure as the even polynomial series. 
Thus, are expressed as 

Too Too TooPr+'^r eXp [-W^+P'i+Pi+Pi] dp^dpydp. 



a: 



odd 



Too Too I-JS°J'y exp [-C^l + Pi + Pi + Pi] dp.Jpydp. 



(B17) 



From this equation and the relation (B6 ), one can derive the expression of the odd polynomial 
series. The first three are 



S, 



odd 



5° 



odd 



Pz, 
Pi- 



^odd 



Pz 



;Pz 



C 5^4(0^2(0 -3[i^3(CP 



Pz 



3i^3(C) 



Pz 



(B18a) 
(B18b) 

Pz- (B18c) 



4. Properties 



Here I list some properties of the two orthogonal polynomial series derived above. Because 
the products S'^^^S^'^ are odd function with respect to Pz, the following integrals vanish; 



oo poo poo 



ceven oodd 

o„ exp 



oo J — oo J — oo 
oo poo poo 

even codd 



CJ'^+Px+Py+Pz 



dp^dpydpz 



oo J —oo J —oo 



ceven oodcl 



-C\li + Pl+Pl + Pl 



dp^dpydpz 



^Ji+pI + pI + pI 



0, (B19a) 
0. (B19b) 



Next, the derivatives of the polynomials can be expanded by 5"° 



odd. 



n-1 



dpz 



ED even O' 



odd 
k ■ 



(B20) 



fe=0 



Because the derivatives of ST'^'^ is an odd function. On the other hand, the derivatives of 



5*° can be expanded by 5*^™°; 



dS: 



odd 



Epodd Q( 
^nk ^/ 



even 
k 



(B21) 



fc=0 



The expressions of -B^X'^" and -8°^'^ are obtained by directly differentiating the expressions 



(B14), (BIS), and their extensions for larger n. 
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FIG. 1: The configuration of the two-dimensional relativistic Vlasov-Maxwell equilibrium derived 
by applying the method in this paper. The gray scale represents the density distribution of electrons 
(or ions) in arbitrary units and the arrows represent the magnetic fields. 
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FIG. 2: (Color online) The momentum distributions of ions at the 0-point (solid line) and the 
X-point (dashed line). 
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